Chiral Limit of Nucleon Lattice Electromagnetic Form Factors 

Walter Wilcox 
Department of Physics, Baylor University, Waco, TX 76798 

Terrence Draper and Keh-Fei Liu 

Department of Physics, University of Kentucky, Lexington, KY 40506 



We calculate electric and magnetic form factors of protons and neutrons in 
quenched Monte Carlo lattice QCD on a 16 3 x 24 lattice at f3 = 6.0 using Wil- 
son fermions. We employ a method which characterizes one of the nucleon fields as 
a fixed zero-momentum secondary source. Extrapolating the overall data set to the 
chiral limit, we find acceptable fits for either dipole or monopole forms and extract 
proton and neutron magnetic moments, the magnitude of which are 10 to 15% low 
compared to experiment. In the extrapolation of the dipole fit of the form factors, we 
find that the dipole to nucleon mass ratio is about 7% low compared to experiment. 
In addition, we obtain positive values of the neutron electric form factor, which, how- 
ever, are poorly represented by a popular phenomenological form at intermediate to 
small k values. A zero- momentum technique for extracting hadron magnetic mo- 
ments is briefly discussed and shown to yield unrealistically small magnetic moment 
values. 
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I. INTRODUCTION 



The techniques of Monte Carlo lattice QCD continue to be developed and applied to nu- 
merous quantities of phenomenological interest. In particular, electromagnetic form factors 
are a useful probe of hadron internal structure. The method used in Ref. 1 characterizes 
one of the meson interpolating fields as a zero-momentum secondary source. Here, we apply 
this technique to the nucleon. This allows us, in the final analysis of the quark propagators, 
to reconstruct a number of operator probes, among which are the various lattice axial and 
vector currents. We report here on our results using the conserved lattice vector current; 
see Ref. 2 for preliminary results using an axial probe. Using this technique with separate 
proton and neutron electric and magnetic sources, we extract all four of the nucleon elec- 
tromagnetic form factors at any desired lattice momentum transfer. (For an independent 
implementation of this technique as applied to the proton, see Ref. 3, where, however, only 
the proton electric form factor is studied.) This investigation for the nucleon is complemen- 
tary to those in Refs. 4 and 5 which fix the vector current, rather than the particle field, as 
the secondary source and analyze the full set of spin 1/2 baryon form factors, but only for 
a single momentum transfer value per quark mass value. See also Ref. 6 for an introduction 
to the formalism and Ref. 7 for some preliminary results of the present investigation. 

The purpose of this study is to begin a comprehensive examination of the mass and 
momentum dependence of nucleon form factors in the quenched approximation. In our 
treatment, we will extrapolate the parameters of the functional forms rather than individual 
form factor values. This is required if contact with various phenomenological forms for 
nucleon form factors is to be made. 

We begin with an introduction to the formalism of lattice sources and correlation func- 
tions. We then present our results for the electric and magnetic form factors and their 
extrapolation to the chiral limit. We close with a summary and some comments about the 
directions of future lattice calculations. 
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II. FORMALISM 



A. Fundamentals 



In order to set the notation, we first give a fairly complete catalog of the definitions 
and symbols used in this paper. Our conventions follow Sakurai ||, although we adopt the 
non-standard representation gamma matrices: 
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We use the four-vector notation 6 M = (b,ib ) , where b is purely real, to ease the transition 
to Euclidean space. Our discrete Minkowski fermion action is given by 



I, J 



(3) 



where I = {x, a, a} and J = {y,j3,b}, which are in the order {space-time, Dirac, color}. 
Flavor sums are understood where appropriate. The matrix Mu[W} is defined by 



MuiW] = 6j,j - [(1 - 7„U (U,(x)) ab b xn 
+ (! + 7m) q/ 3 (pl(x - a M )) a& 5 x , y+atx 



(4) 



Using k = l/2(ma + 4) and (all fields not otherwise specified are lattice fields and 'a' is the 
lattice spacing) 



we find that with UUx) = e LaA A x ) ) this corresponds to the continuum action: 



(5) 



Qcont. 



- f d Xl dx 2 dx 3 dt (tp cont \x) ( 7 „ (a„ - zA™ n *-(x)) + m) ^ con *-(x; 



(6) 
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The assigned (1 ± 7^) structure in Eq. (|j) means that the upper components of the Dirac 
equation propagate in the forward time direction in the static (small k) limit. The conserved 
vector current || from Eqs. (|3]) and (f|) is identified from (A M co>(a:) = uj(x + a M ) — ui(x)) 

$(x) -»• e~ iu3 ^i;(x), (7) 

4>(x) -> ^(x)e iuj{ - x \ (8) 
5Sp 

J »- 5[A,co(x)Y W 
— 0') is given explicitly by (for a single flavor) 

J M (x) = m fax + a„) (1 + 7^) UfoWx) - ^(x) (1 - 7m ) U^x^x + O) . (10) 

p(x) is normalized by (t > t' > 0) 

^(vac|<(x)p(x')^(0)|vac) = <vac|<(x)^(0)|vac>. (11) 

The interpolating fields we use for the proton are 

Xa (x) = e abc ^\x)4 ) \x){C)^{x), (12) 
xAx) = -6 a ' b ' c '^ )c '(C) y/3 ,^ )6 '(x)^ )a '(z), (13) 

where u <-> d for the neutron. (We assume m u = throughout.) The charge conjugation 
matrix C = 72 satisfies C^C -1 = 7* where C = C75. We now continue to Euclidean space 
(t — —it, where £ remains real) and use the integration formula [|H| 



(vac|T(V a (-it A )^(-^) • • Olvac) = Z' 1 J dUd(d(e- s a- sE F ^( a (t A )at B ) . . . , (14) 

where the (, ( are independent, totally ant i- commuting Grassmann integration variables and 
Sq, Sp[(,C] are the Euclidean gluon, fermion actions [[□]]. For example, on our lattice (a 
and (3 are generic indices) 

/ d(d( CaC/3 e~ ZM( = det(M)S aP , (15) 
/ d(d(( a ( p ( 7 ( s e-< M t = det(M)(S aP S 7S - S a5 S, p ), (16) 

where S = M -1 . 



B. Source technique 



For the proton two-point function, we define (assuming t > in the second line) 



G pp (t;p,T) = Y,e- ifS r a , a (™c\T( X a(x)Xa>m I vac) 



= E 



« **r a , a e^(-e° w )(<^(<^ 



I 1 ? 



(vac|V£ u)a (x)4 u)6 (x)^ )c (x)V7 C ' (0)^T (O)^T' (0)|vac). 



(17) 



We will specify later the particular T-matrices which we use here and in following equations. 
Defining in Dirac space the quantity Q = (CQC~ 1 ') for an arbitrary matrix Q, we have ( 
setting det(M) = 1) 



G pp {t;p,T) = Y j e- i ^e abc e a ' h ' c ' (tr [TS (u)aa ' {x,0)S (d)bb ' {x, 0)S {u)cc '{x, 0) 

X 

+tr [TS^ aa '{x, 0)] tr [S {d)bb '{x, 0)S^ cc '{x, 0)]) , 



(18) 



where a configuration average is understood and the trace is over Dirac indices only. The 
corresponding connected Wick contractions of the proton three-point function are given as 
(q = p-p') 

G pJlMP (h, t i; p,p', r) = -i £ e-^^e^T a , a (vac|T (XcM^MxciO)) |vac) 
= Y / e- t v- S2 e abc e a ' b ' c ' (tr [rg ul S (u)aa '(:r 2 , 0; t u q, ^)S {d)W (x 2 , 0)S^ u)cc '(x 2 , 0) 



X2 



+ tr 



TS^ aa '(x 2 , 0)q d s! d)W \x 2 , 0; h, q, fi)S^ cc '(x 2 , 0) 
+ tr [TS^ aa ' (x 2 , 0)S {d)W (x 2 , 0)q u S^ cc ' (x 2 , 0; h, <f, fi) 
+ tr [rg u ^") aa '(a; 2 ,0;ti,g,/i)] tr [S {d)bb ' (x 2 , 0)S^ cc \x 2 , 0) 



+ tr [r5-(«)«»'(a; 2 ,0)]tr g d £ U (x 2 , 0; t u q, fi)S^ cc '(x 2 , 0) 
+ tr [r^ aa '(a; 2 ,0)] tr [S (d)W (x 2 , 0)q u S^ cc ' (x 2 , 0; t u q, //)]) . 



(19) 



This is similar to the two-point function, except that each quark propagator, 5, has been 
replaced, in turn, by qfS where qj is the quark charge ( q u = 2/3, q d = —1/3) and 
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5(z 2 ,0; = K J2e im (S(x 2 ,x 1 + a fl )(l + lfl )Ul(x 1 )S(x 1 ,0) 



- S(x 2 , xi)(l - y^)U^{x 1 )S{xi + a M , 0)) 



(20) 



describes the quark propagator coupled, with momentum q, to the electromagnetic current 
given in Eq. (|T0"D. In order to compactify our notation, let us introduce the quantity 



(SA(x 2 ,o))y a 



e e 



+ tr 



(s bb '(x 2 ,o)s cc '(x 2 ,o)r) , +(rs bb '(x 2 ,o)s cc '(x 2 ,o)) , 

V /a' a \ / a c 



S bb '(x 2 , 0)S cc '(x 2 , 0)1 T a , a + tr \TS bb '(x 2 , 0) 



s c 



(21) 



corresponding to the u (d) quark contribution to the three-point function for the proton 
(neutron), as well as 



(S B (x 2 ,0))Z a a =e^ Vd [(S bb '(x 2 ,0)TS cc '(x 2 ,0)) ai 
+ tr\TS bb '(x 2 ,0)} (S cc '(x 2 ,0)) , 



(22) 



corresponding to the d (u) quark contribution to the proton (neutron). In addition, define 



S(x 2 ,0;tx,q,fi) = Y^S(x 2 ,y)X(y,0;t 1 ,q,n), 
y 

where X(y, 0; t±, q, //) is given explicitly as 

- S V)X1 (1 - y^)U lJL (x 1 )S(xi + 6^,0)]. 
Then we may write concisely for the proton and neutron 



(23) 



(24) 



Gpj^fatDPif \Y) = quAjpfahifip^r) + q d B Jii (t 2 ,t 1 ;p,p',r), (25) 
G nJfl n(t2,ti;p,p',T) = qtAj fahifip^T) + q v Bj fahifip^T), (26) 



where 



Aj^M-JJ'^) = J2e- ifx '' 2 Tv[S A (x 2 ,0)S(x 2 ,y)X(y,0-,ti,q^)}, (27) 
Bj^t 2 M-JJ'^) = Y,e^ a Tr[S B (x2,0)S(x 2 ,y)X{y,0-,ti,q,ri]- (28) 



The 'Tr' notation denotes a trace over both color and Dirac indices. 

We have succeeded in rewriting the three-point functions for the proton and the neutron 



quantities can not be calculated directly because of the presence of the S(x2, y) propagators, 
since there are sums over both x 2 and y present. These equations also make it clear that 
there are two remedies for this situation. One possibility is to introduce a source to simulate 
the current, contained in the X(y, 0; t±, q, /i) factor and associated with the y sum above. 
This technique, of course, is not specific to the nucleon and works for any hadron field. 
However, as one can see from Eq. fl24|) , this choice fixes the spatial momentum transfer, q, 
for a given set of quark propagators. The other possibility is to introduce a source to simulate 
the two quark lines, contained in Sa and S B , which lead to the final nucleon associated with 
the £2 sum. This technique is specific to the nucleon, but leaves the spatial momentum 
transfer free. Indeed, even the choice of which operator to reconstruct is deferred until the 
final analysis, so that the propagators calculated are also useful in studying, for example, 
an axial current. In summary, fixing the current allows one to do a survey of particles with 
a given probe, whereas fixing the particle source allows one to use a variety of probes on a 
single particle. We choose here to fix the particle source since we are particularly interested 
in trying to understand the q 2 dependence of the nucleon form factors. 
For this purpose we introduce v T A B (y]t2,p) such that 



in a very compact 




make it very clear that these 



(29) 



y 



Bj,Xt2,ti;p,p',T) =J2 Tt [ v B(y^2,P)X(y,0;t u q,fi) 



(30) 



y 



where the transpose is over both Dirac and color indices. Explicitly, 




(31) 



We now multiply on the right of Eq. (pHf ) by M(y, x') and sum on y to give 



J2vl B (y;t 2 ,p)M(y,x') = e'^' S A , B (x', 0) . 



(32) 



y 
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Using the well-known relation 

M\x,y)= l5 M{y,x)j 5 , (33) 

where 'f works in Dirac and color space, one can show that Eq. fl32|) leads to 

^M(x',y) l5 vl B (y;t 2 ,p) = e^Vi.B^', 0) . (34) 
y 

The right hand side of fl34|) identifies the sources 

fe^B^; t 2 ,p) = ^ s ' l5 S\ B {x', 0) , (35) 

to be inserted in the matrix inverter. Actually, for our purposes, it is more convenient to 
consider a linear combinations of sources, 

W roton = q u b A + q d b Bl (36) 
b neutron = q d b A + q u b B , (37) 

which give the desired proton and neutron three-point functions directly. (See Eqs. ( fZ5D and 
above.) 

The general method of using secondary or sequential sources to perform spatial sums 
over intermediate lattice operators was introduced in Refs. 12. 

C. Correlation Functions 

Although the behavior of Euclidean-time nucleon correlation functions on the lattice is 
standard material 0-0, we summarize their properties for future reference below. 

We deal with both protons and neutrons, but for simplicity the following discussion will 
be for the proton only. The two-point function defined in Eq. ( |T7D in the large Euclidean 
time limit gives 

G pp (t-p,T) 7V^ e - B *r a , a (vac|xa(0)|^s)(^s|x a ,(0)|vac), (38) 
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where N v is the number of spatial points in the lattice volume. We write lattice and 
continuum completeness, respectively, for fermionic states as 



J2\n,P,s){n,p,s\ = I, 



n,p,s 



The continuum limit 



V- / d 3 p m , 



d 3 p 



(39) 
(40) 



V 



E 



(2tt)* 



(41) 



where V = N v a 3 , gives the correspondence between lattice and continuum states: 



\n,p,s) — 
Use of Eqs. (|5|) and ( f42"|) then give 
(vac\xa(0)\p, s) - 

(P>|Xa'(°)l VaC ) - 



m \ 1 I 2 

ve) lw 



(42) 



m 



1/2 



(2k) 3 / 2 ViV„£, 

3 / „ x l/2 



vac|xr"(0)b>), 



m 



(P,s|r^(0)|vac). 



(43) 
(44) 



(2k) 3 / 2 ViV„E, 

(vac|Xa° nt '(0)|p, s) and (p, s|%y n *'(0)|vac) are related as usual to the free spinors u a (p, s) and 
u a '(p,s) by 



(vac|x™ n *'(0)|p, s) = Z?i a (p, s), 
(p,s\T°> nt iO)\™c) = Z*u a ,{p,s), 



(45) 
(46) 



where the unknown amplitude, Z, transforms as a scalar and therefore can not depend upon 
p or s. Now, using these results in Eq. (|38|) , one obtains 



G pp (t;p,T) — > - — — e tr 



-i ® + m 



(2k) 3 £ 

where we have used the following relation for free spinors: 

[— i p 1 + m) 



2m 



(47) 



Ys u iv, s)u(p, s) 



2m 



(48) 
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For example, with the choice 



r 4 = 



/ 




one obtains 



2E (2k) 3 

For the proton three-point function, Eq. ([19|) , the large time limit is 



(*2-ti),ti>l ;v2\- ,-E p {t2-H) e -E p ,ti 

r Q '„(vac|xa(0)|p,s)(p, S |J M (0)|p',s'}(p / ,s / |Xa'(0)|vac). 



The lattice and continuum matrix elements of J M are related by 

,2 \l/2 



<£s|J M (0)|jf,s') 



m 



(p,s\r° nt -(o)\p',s'), 



where ( a^ u = hilmlv] and F\ and F 2 are real) 



(p, ^1^(0)1^, s') = iu(p, s) (7^ - <v|^*a) M G°~"> 

Based on these forms, with the choices p = 0, /x = 4 and T = T 4 , the proton three- 
function yields 

<W*». H 0, -« r 4 ) ^ »' l|^ e -*-*) e -«. (^2) c.fo"), 

where G e (q 2 ) = [F\ — ( 2 m) 2 ^ 2 ) * s ^ e e l ec tric form factor. Similarly, using 



r, = - 



(Tfe 





one finds that with p = 0, /i = j and T = 

2£ (2k) 3 

where G m = (Fi + F 2 ) is the magnetic form factor 
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Eqs. (|50|), (p4|), (|5q) provide the means of calculating G e and G m . We do not measure 
these correlation functions directly, but instead analyze the ratio 

r(*)= ( 2E V"^ f G ^W2;Q,-g,r4) \ ( G pp (t'-,o,r 4 ) \ & 

e[q ] \E + m) \ G pp (t 2 ;0,r 4 ) ) \G pp (t>; q, T 4 ) ) ' 1 j 

for G P and 



1 2^ ^ + m \ 2F ( G pJjP (t 2 ,t 2 /2;0,-q,T k ) \ ^ G pp (t'; 0, T 4 ) 
' " 2 £j % I 2£ J 1, G pp (t 2 ;0,r 4 ) A ^5 £ I*) 1 ' ^ 



'2 
2(' 



for G m , which are seen to have smaller error bars. Notice that Eq. (p7j) guarantees that we 
measure G e (0) = 1 except for convergence errors and that all currents are located halfway 
between the creation and annihilation time steps. Also notice that the fractional powers 
are present in order to provide a choice in the two point functions used in order to allow 
single exponential behavior to develop in t'. The identities developed in Refs. 1 and 6 show 
that Gp P (t; q, T 4 ), G p j 4P (t 2 , t\\ 0, — q, T 4 ) and G p j jP (t 2 , t\, 0, — q, T^) are real functions in the 
ensemble average after evenness (Eqs. ( |5T)D and and (|M[)) or oddness (Eq. flSTSP) in spatial 
momentum <f has been enforced by hand. Because the time interval between the initial and 
final nucleon sources here is odd, t 2 = 15, we define a spatial current operator at half-time 
steps according to 

Jj(xi,t x + -) = - {Jj{x u t x ) + Jj{x u h + 1)) (59) 

for use in (|56D . The time-nonlocal charge density operator, J 4 , is already naturally associated 
with half-time steps. 

III. RESULTS 

Our quenched configurations are 16 3 x 24 and were calculated using the Monte Carlo 
Cabibbo-Marinari pseudoheatbath [HJ. The SU(3) fundamental Wilson action was used 



with periodic boundary conditions and (3 = 6.0. The gauge field was thermalized for 5000 
sweeps from a cold start and 12 configurations separated by at least 1000 sweeps were saved. 
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Since our gauge configurations are taken from a single Markov chain, if the resulting correla- 
tions between configurations are sufficiently small then the configurations chosen for analysis 
are effectively statistically independent. The preferred way to test for autocorrelations is to 
bin consecutive configurations; we have too few configurations (12) for this to be a useful 



test. A binning of a larger data set [I4~ on quenched lattice at (3 = 6.2 revealed no evidence 
of autocorrelations for two-point hadronic functions when configurations were separated by 
250-500 (multi-hit Metropolis) sweeps. An examination of a number of observables as a 
function of Monte Carlo sweep number in the quenched (3 = 5.9 form-factor calculation of 
Ref. 5 revealed no apparent auto-correlations for a sample of 28 configurations separated by 
1000 pseudo-heatbath sweeps. The same sweep-number separation was chosen for calcula- 
tions of semi-leptonic form factors in Ref. 15; we have made this same conservative choice 
in the present calculation to keep auto-correlations to acceptably low levels. 

For the quarks we use periodic boundary conditions in the spatial directions and 'fixed' 
time boundary conditions, which consist of setting the quark couplings across the time 
edge to zero. The origin of all quark propagators was chosen to be at lattice time site 5; 
the secondary zero- momentum nucleon source was fixed at time site 20. We expect that 
these positions are sufficiently far from the lattice time boundaries to avoid nonvacuum 
contaminations. All our results for proton form factors use the point interpolating fields, 
Eqs. fljjD and ([13]), and similarly for the neutron. We used the conditioned conjugate gradient 
technique for quark propagator evaluation described in Ref. 16. For our convergence criterion 
we demanded that the relative change in the absolute sum of the squares of the quark or 
secondary propagators be less than 5 x 10~ 5 over 5 iterations. As one check of the nucleon 
secondary source, we verified current conservation for t2 > t\ > to (9(10~ 4 ). Since we 
wish to calculate the electric and magnetic form factors of both the proton and the neutron, 
one non-source and four source propagator inversions are necessary per configuration. The 
results below include k = .154, .152, .148, and .140. Our statistical error bars come from 
first (form factors), second (form factor fits) and third order (chiral extrapolations) single- 
elimination jackknifes. 
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We show the pion, rho and nucleoli masses measured on our 12 configurations in Table I. 
(The interpolation fields used for the pion and rho are the usual relativistic ones.) Actually, 
we show the results of single exponential mass fits using both smeared (over the entire 
lattice volume using the lattice Coulomb gauge) and point quark propagators. The fits here 
are for lattice time sites 16 to 19 for k = .154,. 152 and .148 and time sites 18 to 21 for 
k = .140, which was seen to evolve more slowly in time. The smeared and point masses 
are consistent with one another within the statistical error bars, but a small systematic 
downward shift of the smeared masses relative to the point masses seems to be present. 
The smeared results are also consistent with the more accurate (3 = 6.0 mass results of 
Ref. 17, with which our largest two k values overlap. When needed, the smeared mass 
results from Table I (from our simulation) will be used for the kinematic factors which 
appear in (^) and (|58"D; the uncertainties associated the kinematics are then included in the 
form factors as uncorrelated errors, which affects mainly the magnetic error bars. We will 
also use the accurately determined (3 = 6.0 value of n cr = .15708(2) from Ref. 18 for our 
chiral extrapolations. In the following, we will illustrate our data mainly with the k = .154 
results, where, it must be kept in mind, the error bars are the worst. 

In order to test for continuum dispersion and to examine the time dependence of our 
two point functions, we define the local mass, energy, and energy minus mass from (|50|) as 



These quantities are given in Fig. 1 for k = .154 as a function of lattice time. The starting 
position of all quark propagators is time step 5. The horizontal lines in this figure give the 
expected result from the continuum dispersion relation using the central value of the mea- 
sured nucleon mass. The most significant quantity relevant to our form factor measurements 




(60) 




(61) 




(62) 
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is Eq. (|62|) since this involves a ratio of quantities that enters in Eqs. fl57|) and ([58]). From 
the figure, it appears that this measurement does not become consistent with continuum 
dispersion until time slice 14|, which corresponds to t' — 9 in Eqs. fl5"T|) and (|58"D. (The 
propagator time origin is defined to be t' — 0). The other k values are similar, and t' — 9 is 
used in all of our form factor results. 

Figures 2 and 3 represent energy minus mass measurements as measured from ratios of 
the three-point functions, Eqs. (|54]) and (|56|) . These are defined at integer time steps (the 
currents at t\ are defined at half-integer sites) by 

- ml(tl + 5>l* = ln U*»ft,«. + i;0,-?,rj • (63) 

G p j k p{t 2 ,ti]Q, -q, Tj) \ 
Gpj^h + l-O,-^)) • 1 ; 

Figure 2 shows the local [E—m] values from the magnetic three-point function. This function 
is seen to have a flat exponential behavior which begins quite near the source origin. Figure 3 
shows the [E — m] values from the electric three-point function, which, unlike the magnetic 
case, never appears to flatten. On the other hand, the results are consistent with the 
expected [E — m] values near the mid-point of the lattice where the data for Eqs. ( |57| ) and 
(|58l) are actually taken. 

The results of our form factor measurements are presented in Tables |lj through 0. 
Figures 4 though 7 present the k = .154 graphical results. Notice that the (qa) 2 > 
values given in the tables are affected very little by the error bars in the nucleon mass (the 
maximum error bar, at k = .154, is approximately 3% of the central value). Our philosophy 
in comparing our results to experiment is to look at the simplest phenomenological forms 
consistent with the lattice data, and then to extrapolate the fit parameters, rather than the 
individual form factor values, to the chiral limit. The solid and broken lines in these graphs 
represent the best simultaneous dipole and monopole fits, respectively, of the combined 
proton electric, magnetic and neutron magnetic form factors. These are three parameter 
fits, giving the fit dipole mass from 
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[E-m](h + -)\ jk =]n 



G D M) = - ^r, (65) 

1 + 



i 2 



D 



or the monopole mass from 

G? (q 2 ) = —V > ( 66 ) 
as well as the proton and neutron magnetic moments from the forms, 

G D m {q 2 )= Gm ^\ 2 , (67) 
1 + 



or 



m D 



G A M) = . (68) 

The fit parameters found this way are listed in Table [V| along with the chi-squared per 
degree of freedom found (Xd)i and the dipole/nucleon or monopole/nucleon mass ratio. In 
general, the quality of the fits are reasonable except for the monopole form at k — .140. 
The simultaneous monopole fits are seen to be slightly preferable to the dipole ones at the 
lowest three k values. Notice also that the ratio tum/itin is essentially flat over these three 
k values. 



In a separate fit of the proton electric form factors, we list in Table [VII| the dipole and 
monopole masses, the corresponding charge radii (in lattice units), mass ratios, and x\ values 



found. In comparing the x\ values from the dipole fits in Tables [\T| and |VII| , we notice that 



the inclusion of the magnetic data decreases the quality of the fits at k — .154 and .140, 
whereas it increases the quality at the two intermediate k values. For the monopole fits, we 
see that the x\ values are low in all cases, except again for the k = .140 simultaneous fit. It 
is the inclusion of the magnetic data that is responsible for the large x\ value; the values in 
Table |V| inform us that the (qa) 2 falloff of the magnetic data is faster than for the proton 
electric data at this k value. 

Figures 7 through 10 show the measurements of the neutron electric form factor at the 
four k values. The phenomenological form, 
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G n M 2 ) = ~^rG n m {q 2 ) , (69) 

is compared to the numerical results, using either the dipole (solid line) or monopole (broken 
line) parameters from Table [VTl to characterize G^(g 2 ). Although the data are quite noisy, 
we obtain positive values of G™(q 2 ) in all cases, in agreement with most experiments in 
this energy regime. In addition, we see that the above phenomenological form fails at the 
two lowest k values. This does not rule out such a form in the chiral limit, but does make 
whatever physics lies behind it less compelling. 

Figure 11 shows the k = .152 proton electric form factor. In our measurements we have 
the option of reconstructing the spatial momentum transfers in a number of different direc- 
tions for a given (qa) 2 value. This figure shows the effect on the error bars of averaging (O) 
and and not averaging (□) over equivalent momenta. (The (qa) 2 values of the nonaveraged 
data have been increased slightly so that the two data sets do not overlap.) At these values 
of k the effect is to reduce the error bars by a factor of two to three. 

Also shown in Fig. 11 are the best dipole (solid line) and monopole (broken line) fits to 
this data from Table |VI]] . These fits illustrate that the source of the large x\ value seen in 
the dipole fit of Table |VII] at k = .152 (and similarly at .148) comes about because of the 
failure of the highest (qa) 2 measurement to fall off sufficiently fast. It is possible this is a 
systematic high (qa) 2 lattice artifact; on the other hand, the fact that the k = .154 and .140 
results do not display similar behavior undercuts this explanation. 

Figure 12 represents a comparison of the results of two methods of extracting the proton 
magnetic moment. In this figure, the data points given by the square symbols are taken from 
the dipole fits of Table [Vl|, while the diamond data points are given by the zero momentum 
measurement (£2 > £i > 0): 

tr [r fe (vac|x(x 2 )(xi)i J j (xi)x(0)|vac)] 
$2fh _ G m (0) 

]>>[r 4 (vac| X (* 2 )x(0)|vac>] ~ m 2m ' 

The continuum formula on which the above is based is derived by taking the derivative 
of the continuum analog of G p j jP (t2, ti; 0, — q, r&) with respect to the zth component of q, 
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evaluated at q = 0, and dividing by a zero momentum two point function. The resulting 
equation is then transcribed into lattice language by changing the continuum matrix elements 
into lattice ones and by making the substitutions / d 3 x — > aJ2x an d Jf nt \x) — ► a~ 3 Jj(x). 
Although the two results agree at the smallest k, the zero momentum measurements give 
unrealistically small magnetic moments at the larger k values. The reason for this behavior 
is the same as for similar behavior seen in lattice mesons using charge overlap techniques 
|i9| . Because the lattice matrix elements do not contain arbitrarily small momentum states, 
the continuum derivative at q = can not be duplicated, and the lattice version would only 
be expected to hold for D/2 ^> R, where D is the length of the lattice on one side and R 
is a hadron correlation length, say, the charge radius. That is, the hadron is expected to be 
well contained in the given lattice volume. Apparently, this condition is only beginning to 
be satisfied at the smallest measured k value, which is the farthest from the chiral limit. 

Figures 13 and 14 represent the chiral extrapolation of the dipole fit proton and neutron 
magnetic moments. The values found from these fits as well as from similar monopole fits 
are listed in Table |V111| . The quantity m q a is defined to be 

mq<X = 1(k~t)' (71) 

These extrapolations are simply linear fits, which were adequate to describe the data, as 
seen from the Xd values. The magnitudes of the magnetic moments are 10 to 15% low; the 
magnitudes found in Ref. 5 are also low, but by 15 to 30%. There is a hint in these figures 
that the largest k values prefer to lie above the linear fit, so more satisfactory magnetic 
moments may result from a deeper exploration of the region near K, C riUcai where, however, 
error bars more problematical. 

Figure 15 presents the chiral extrapolation of the dipole to nucleon mass ratio from 
Table |V111| . This is assumed to be linear as a function of m q a; again, the \\ values in Tables 



|V111| and |X| do not demand a more sophisticated treatment. We prefer to do the chiral 
extrapolations on the mass ratios from the above monopole and dipole fits because the ratio 
of similar physical quantities is often less subject to systematic errors and because measuring 
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the nucleon mass in relation to other hadrons is best done in a separate high statistics 
spectrum calculation, such as the (3 = 6.0 calculations of Refs. 17 and 18. For comparison, 
we have provided the results of three parameter (Table |V111| ) and one parameter (Table [TXD 
fits of experimental nucleon data taken from Refs. 20 and 21. Setting the scale from the 
chiral limit nucleon mass of Ref. 18, our four-momentum transfer range roughly corresponds 
to .6 GeV 2 < q 2 < 1.9 GeV 2 ; the values overlapping with this range from these references 
are the values used in the fits, which are listed in Tables [X] and XI. The Xd values for the 



experimental results, which are more precise than the lattice measurements, show that more 
parameters are really required in this energy range to produce reasonable phenomenological 
fits. We have \\ = 1.99 and 10.7 for the dipole and monopole fits, respectively, of the 
proton electric data in Table |IX|. Similarly, we obtain \\ = 3-39 and 35.5 for the dipole and 
monopole fits, respectively, of Table [VII]] . In comparing the two experimental fits, we notice 
that the dipole mass ratio is slightly larger when the magnetic data is included (Table |VIII| ) 
than when it is not (Table [IX]). Experiment shows that the proton and neutron magnetic 
form factors fall off significantly more slowly than the proton electric form factor in our 



energy regime [p2| . This explains the tendency of the experimental dipole fit which includes 
the magnetic data to produce a larger value of m^/m^ than a similar fit of the proton 
electric data alone. It is encouraging that the same tendency seems to be present in the 



lattice dipole results in Tables |VIII| and [TXJ, although the overall value for the ratio seems 
to be about 7% low in either case. However, as in the magnetic moment case, the larger k 
values of this ratio prefer to lie above the linear fit, and therefore a more satisfactory value 
of this ratio could result from a deeper exploration of the chiral limit. 

IV. COMPARISON AND CONCLUSIONS 

We have investigated the functional forms of the nucleon electromagnetic form factors 
as given by quenched lattice QCD. Although our results are not sufficiently accurate to 
distinguish between monopole and dipole fits to the data, we have seen that the error bars 
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on these quantities are encouraging and that the magnetic moments as well as the dipole to 
nucleon mass ratio are reproduced to within about 15%, similar to spectrum calculations. 
We have also seen an indication that the chiral limit proton and neutron magnetic form 
factors have a slower falloff in q 2 than the proton electric form factor, which is similar to 
experiment in this energy regime. In addition, the neutron electric form factors come out to 
be positive, but their values are poorly represented by a popular phenomenological form at 
intermediate to small k values. Finally, we have investigated a zero-momentum technique 
for extracting magnetic moments, but found that this method yields unrealistically small 
values approached. 

In agreement with the results of Ref. 5, which used 12 2 x 24 2 lattices at /3 — 5.9, we 
have found that the magnitudes of the proton and neutron magnetic moments are small 
compared to experiment, and, indeed, that the neutron value is more badly represented 
than the proton. The small improvement over the Ref. 5 values could be due to the larger 
(3 value used here. Alternatively, we have seen a hint in Figs. 13 and 14 that perhaps better 
values simply await a closer approach to n cr rather than requiring larger (3 values. Another, 
more interesting, possibility is the unknown contribution to the magnetic moments from 
disconnected quark loops due to current self-contractions The Ref. 5 calculation also 
gave positive values for the neutron electric form factor as was found here. There does 
seem to be a difference in the physical size of the nucleons in these two studies, however. 
Corresponding to the dipole and monopole chiral extrapolations in Table |X|, we find that 
that the dimensionless quantity R p m^ has a value 4.23(18) from the dipole fit and a value 
4.80(35) from the monopole fit. If we simply divide these results by the experimental (aver- 
age) nucleon mass, equivalent to the procedure adopted in Ref. 5, we then obtain a charge 
radius of .89(4)fm from the dipole fit and 1.01(7)fm from the monopole fit as compared to 
the estimate .65(8)fm from Ref. 5. The experimental result is .862(12)fm ||24j| . The Ref. 5 
result corresponds to a dipole to nucleon mass ratio which is larger than the experimental 
result by about 25%, whereas the result here is about 7% low compared to experiment. It is 
unlikely that the physical box size plays an important role in this difference since, assuming 
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renormalization group scaling, the box dimensions of these two calculations are very close 
(comparing with the shortest box dimension of Ref. 5). The systematics associated with 
the different techniques of extracting the form factors could be responsible for these rather 
different results. 

Recent studies of scaling show that both the dimensionless ratios of the string tension K5 



and the scalar glueball mass to the chiral condensate [26[ have scaling violations of ~ 20% 
from ft = 5.7 to ft — 5.9 and ~ 10% from ft — 5.9 to ft — 6.0. Although glueball mass 
studies on large lattices seem to show scaling from ft = 5.9 to ft = 6.2 |27|], hadron masses 
and /„- ratios still show a deviation of the order of ~ 10 to 20% lf28fl . Therefore, masses and 
magnetic moments measured at ft = 6.0 are subject to a scale-breaking systematic error 
which could be as large as ~ 20%, although our use of mass ratios in the monopole and 
dipole fits would be expected to significantly reduce the systematic error in the extrapolation 
of the functional forms to the chiral limit. 

The overall message of form factor measurements to this point seems to be that the 
quenched approximation adequately represents the bulk of the physics of these quantities; 
however, we are still far from being able to test QCD in a precise experimental way on 
the lattice. At the same time, we should keep in mind that another major theme of such 
calculations is the increased physics understanding that will be afforded through increasingly 
sophisticated parametrizations of lattice laboratory data. This has the potential of teaching 
us about the dynamics associated with quark masses, current self-contraction graphs, and 
the quenched approximation. It is clear that significantly larger computer resources will be 
necessary to make substantial progress in our understanding of these issues; we therefore 
look forward to the benefits of improvements in computer technology, such as proposed in 
Ref. 29. 
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TABLES 



TABLE I. Smeared and point propagator mass fits. 



Particle 



k= .154 



.152 



.148 



.140 



pion 
rho 
nucleon 

pion 
rho 
nucleon 

pion 
rho 
nucleon 



.369(9) 
.46(2) 
.74(4) 

.38(2) 
.46(2) 
.73(5) 

.361(1) 
.463(3) 
.721(7) 



Smeared case 
.488(7) 

.54(1) 

.85(2) 
Point case 

.49(1) 

.55(1) 

.87(3) 
Ref. 17 
.474(1) 
.545(2) 
.861(5) 



.677(5) 
.711(9) 
1.12(1) 

.684(9) 
.721(8) 
1.15(2) 



1.015(5) 
1.032(5) 
1.62(1) 

1.027(6) 
1.040(7) 
1.64(1) 



TABLE II. Proton and neutron form factors at k = .154 as a function of lattice four momen- 



tum-transfer squared. 

Form Factor (qaf = .145 .274 .392 .502 

Gf .481(34) .335(34) .223(34) .186(59) 

G p m 1.17(11) .895(90) .69(14) .86(24) 

G n e .033(17) .026(17) .046(38) .068(67) 

G™ -.748(69) -.578(89) -.526(88) -.54(20) 
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TABLE III. Proton and neutron form factors at k = .152 as a function of lattice four momen- 
tum-transfer squared. 



Form Factor (qa) 2 = .147 .281 .406 .522 

Gl .551(18) .391(24) .293(29) .288(32) 

G p m 1.22(7) .906(59) .696(95) .65(13) 

G n e .0230(77) .019(11) .026(20) .038(29) 

G™ -.781(59) -.586(47) -.474(49) -.406(90) 



TABLE IV. Proton and neutron form factors at n = .148 as a function of lattice four momen- 
tum-transfer squared. 

Form Factor (qa) 2 = .150 .291 .426 .555 

Gl .665(10) .496(17) .390(25) .367(27) 

1.42(4) 1.03(5) .795(66) .643(77) 

G n e .0109(31) .0116(57) .0129(81) .016(10) 

G™ -.907(36) -.666(32) -.520(39) -.403(49) 



TABLE V. Proton and neutron form factors at k = .140 as a function of lattice four- momentum 
transfer squared. 

Form Factor (qa) 2 = .152 .300 .444 .584 

Gl .791(8) .646(14) .539(18) .484(21) 

1.44(3) 1.12(3) .895(41) .740(37) 

G n e .0030(10) .0043(16) .0048(22) .041(28) 

G™ -.927(22) -.726(25) -.582(30) -.476(22) 
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TABLE VI. Combined dipole and monopole fits of the proton electric and magnetic, neutron 
magnetic form factors. 







re — 1 54 


.152 


148 


140 


Dipole 






.688(21) 


noo/ i o\ 






II I'D 1 1' IN 




80Qf31 ) 


74^ f 17^1 








9 ^1 (, 


9 99f1 0) 


9 0Qf6) 








— 1 59143) 


— 1 4514 0) 

-L . 1 1UI 


— 1 34f5) 






Xd 


1.36 


1.45 


.66 


1.63 


lYTononolp 


777, A/f n, 


364f31) 


.427(19) 


539(16) 


725f24) 




m M /m N 


.492(49) 


.503(26) 


.481(15) 


.447(15) 




G p m (0) 


2.63(22) 


2.26(10) 


2.09(6) 


1.77(6) 




G™(0) 


1 R8f 1 S^l 
— l.Do^lo J 




— 1 .OO^O ) 






xl 


.70 


.25 


.64 


3.57 


TABLE VII. Dipole 


and monopole fits 


of the proton 


electric form factor. 








re = .154 


1 

.loz 


.14o 


.141) 


Dipole 




593f32) 


679f20) 


839(4 6) 


1 12CV) 




Rp/a 


5 84f32) 


5 inn 5"i 


4 1 6f8) 

^.1UI (J J 


3 09(8") 




m D /m N 


8D1 f61) 


799f30) 


743(16") 

• 1 TT'.J 1 1UI 


688(19) 




Xd 


.25 


2.24 


1.97 


.41 


Monopole 


mMO, 


.359(31) 


.426(18) 


.542(16) 


.744(23) 




Rp/a 


6.82(59) 


5.75(24) 


4.52(13) 


3.29(10) 




m M /m N 


.485(49) 


.501(25) 


.484(15) 


.459(15) 




xl 


.30 


.41 


.42 


.49 
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TABLE VIII. Chiral extrapolation of the proton, neutron magnetic moments and the mass 









value 


x\ 


Exp'l fit 






ran Itun 


.836(31) 


.39 


.900(5) 


Dipole 




G p m (0) 


2.44(11) 


.64 


2.87(4) 






G™(0) 


-1.59(11) 


.23 


-2.02(5) 






m M /m N 


.516(31) 


.06 


A A O / A \ 

.448(4) 


Monopole 




G p m (0) 


2.48(10) 


.81 


3.22(6) 






G£(0) 


-1.59(12) 


.51 


-2.33(6) 


TABLE IX. 


Chiral extrapolation of the dipole and monopole proton mass 


ratios from Table VII. 








value 


x 2 d 


Exp'l fit 






m D /m N 


.818(35) 


.31 


.883(6) 






m M /m N 


.510(37) 


.08 


.472(4) 



TABLE X. Experimental data for the proton from Ref. 20. 



g 2 [(GeV/cf] G^f) G m (g 2 ) 

.65 .265(12) .767(8) 

.72 .270(17) .690(10) 

.78 .217(7) .647(6) 

.94 .196(8) .523(6) 

1.1 .141(5) .452(4) 

1.35 .114(5) .352(3) 

1.75 .0713(64) .248(3) 
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TABLE XI. Experimental data for the neutron from Ref. 21. 



q 2 [(GeV/c) 2 } G m (q 2 



.60 -.629(20) 

.78 -.434(23) 

1.0 -.345(27) 

1.0 -.322(18) 

1.17 -.284(26) 

1.53 -.203(11) 

1.80 -.185(17) 
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Figure Captions 

1. Local energy, mass or energy minus mass for the nucleon two point function, Eqs. fl60|) 
through (|62| ) for k = .154 as a function of lattice time slice. In this and the next 
two figures, the horizontal lines give the expected continuum results based upon the 
measured mass. All propagators begin at time slice 5. 

2. Local energy minus mass for the proton magnetic three point function, Eq. ( |64|) at 
k = .154 as a function of lattice time slice. 

3. Local energy minus mass for the proton electric three point function, Eq. (|63D at 
hi = .154 as a function of lattice time slice. 

4. Electric form factor of the proton at k = .154 as a function of lattice Minkowski four 
momentum transfer squared. The solid and broken lines represent the best simultane- 
ous dipole and monopole fits, respectively, of the proton electric, magnetic and neutron 
magnetic form factors. The parameters of the fits are listed in Table [VT|. 

5. Magnetic form factor of the proton at k = .154 as a function of lattice Minkowski four 
momentum transfer squared. The solid and broken lines are simultaneous dipole and 
monopole fits, respectively, as explained in the Fig. 4 caption. 

6. Magnetic form factor of the neutron at k — .154 as a function of lattice Minkowski 
four momentum transfer squared. The solid and broken lines are simultaneous dipole 
and monopole fits, respectively, as explained in the Fig. 4 caption. 

7. Electric form factor of the neutron at k = .154 as a function of lattice Minkowski four 
momentum transfer squared. The lines shown represent the phenomenological form, 
Eq. (|69|), using either dipole (solid line) or monopole (broken line) fits of the neutron 
magnetic form factor from Table (VJ. 

8. Same as Fig. 7, except at k = .152. 
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9. Same as Fig. 7, except at k = .148. 

10. Same as Fig. 7, except at k — .140. 

11. A comparison at k = .152 of the proton electric form factors calculated by averaging 
all equivalent momentum transfers (O) and by not averaging (□). The solid and 
broken lines represent the best dipole and monopole fits of the proton electric factor, 



respectively, from Table [VTl 



12. The Magnetic moment of the proton at four values of m q , which is defined by Eq. (fTT]). 
The values from the simultaneous fits given in Table [VI| (□) are contrasted with the 
extracted values from Eq. (|70|) (O). 



13. Linear chiral extrapolation in m q of the proton magnetic moment from the simultane- 
ous fits of Table [VT]. 

14. Linear chiral extrapolation in m q of the neutron magnetic moment from the simulta- 
neous fits of Table [VI]. 

15. Linear chiral extrapolation in m q of the dipole to nucleon mass ratio from the simul- 
taneous fits of Table [Vj. 
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